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ABSTRACT 

One of the most important questions regarding the progenitor systems of Type la supernovae (SNe la) is 
whether mergers of two white dwarfs can lead to explosions that reproduce observations of normal events. 
Here we present a fully three-dimensional simulation of a violent merger of two carbon-oxygen white dwarfs 
with masses of O.9M and 1.1M© combining very high resolution and exact initial conditions. A well-tested 
combination of codes is used to study the system. We start with the dynamical inspiral phase and follow the 
subsequent thermonuclear explosion under the plausible assumption that a detonation forms in the process 
of merging. We then perform detailed nucleosynthesis calculations and radiative transfer simulations to pre- 
dict synthetic observables from the homologously expanding supernova ejecta. We find that synthetic color 
lightcurves of our merger, which produces about O.62M of 56 Ni, show good agreement with those observed 
for normal SNe la in all wave bands from U to K. Line velocities in synthetic spectra around maximum light 
also agree well with observations. We conclude, that violent mergers of massive white dwarfs can closely re- 
semble normal SNe la. Therefore, depending on the number of such massive systems available these mergers 
may contribute at least a small fraction to the observed population of normal SNe la. 
Subject headings: stars: supernovae: general - hydrodynamics - binaries: close - radiative transfer 
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1. INTRODUCTION 

Despite many years of dedicated research, the progenitor 
systems and explosion mechanisms of SNe la remain un- 
clear. It is generally accepted that SNe la originate from 
thermonuclear explosions of massive carbon-oxygen white 
dwarfs in binary systems. Depending on the nature of the 
companion star, two different progenitor sys tems have been 
propo sed. In the single degenerate scenario ( Whelan & Ibenl 
119731) a carbon-oxygen white dwarf accretes from a non- 
degenerate companion star until it reaches the Chandrasekhar- 
mass and explodes (but note that an explosion befo re reaching 
the C h andra sekhar mass may also be possible, e.g. Fin k et all 
I2007L 120 Id). In con trast, in the double degenerate sce- 
nario (Iben & Tutukov 1984; Webbink 1984) a merger of two 
carbon-oxygen white dwarfs causes a thermonuclear explo- 
sion of the merged system. 

Many recent fi ndings including S N rates from population 
synthesis studies dRuiter et al .120091) . s tudies of the delay time 
distribution of observed SNe la (e.g. iMaoz et all 120101), th e 
lack of radio emission of SNe la (e.g. lHancock et al .11201 ll) . 
the lack of hyd rogen emission in nebular spectra of SNe la 
([Leonard 2007 ), the lack of X -ray emission in elliptical galax- 
ies (Gilfanov & Bogdan 2010) and studies of SN la remnants 



(e.g. lBadenes et al.ll2007|) seem to favor the double degener- 
ate scenario. In addition, first studies of SN 20 1 lfe seem to 
disfav or a single degenerate progenitor (see, e.g. lBloom et al.l 
1201 ll) . There is, however, no unambiguous proof for any 
progenitor scenario yet. For a detailed discussion abo ut con- 
straints on the progenitor scenarios see lHoweil (1201 1|) . 

To date, most theoretical work on SNe la has concentrated 
on the single degenerate scenario, as mergers of white dwarfs 
were thought not to lead to thermonuclear explosions. This 
was mainly based on the picture that mergers would leave be- 
hind the more massive white dwarf with a hot envelope made 
up of the material of the less massive white dwarf. Fast accre- 
tion from the envelope will then turn the carb on-oxygen white 
dwarf into an oxygen-neo n white dwarf dNomoto & Ibenl 
ll985l:ISaio & Nomotol 19981) which collaps es to a neutron star 
as it a pproaches the Chandrasekhar mass dNomoto & Kondol 
1199 1|) . Ways to avoid the transformation of the carbon-oxygen 
white dwa rf into an oxygen- neon white dwarf have been pro- 
posed (e.g. lYoon et al.ll2007|) . but without conclusive results. 

Recently, however, we demonstrated that violent mergers 
of two carbon-oxygen white dwarfs could directly lead to 
a thermonuclear exp losion while the merger is still ongoing 
dPakmor et al.ll2010|) . We also showed that the observables 
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Fig. 1. — Snapshots of the merger of a 1.1 M© and a 0.9 Mq carbon-oxygen white dwarf and the subsequent thermonuclear explosion. At the start of the 
simulation the binary system has an orbital period of « 35 s. The black cross indicates the position where the detonation is ignited. The black line shows 
the position of the detonation front. Color-coded is the logarithm of the density. The last two panels have a different color scale ranging from 10 _4 g cm -3 to 
10 6 g cm -3 and 10 4 g cm -3 , respectively. 

for such an explosion with two white dwarfs of 0.9 M show 
good agreement with the observed properties of subluminous 
1991bg-like SNe la. Furthermore, we found that for a primary 
mass of 0.9 M a mass ratio of at least abou t 0.8 is required 
for the sc enario to work dPakmor et aT1l201 ll) . 

Lately, iDan et al.l (1201 ll) showed that using exact initial 
conditions can change the properties of the merger. In partic- 
ular, this leads to a significantly lon ger inspiral in their simu- 
lations. However. IDan et all (1201 1|) were only able to run the 
simulation with a resolution of 2 x 10 5 particles (for compar- 
ison, the violent merger calculations by Pak mor et al.l (1201 Oi 
12011b used 2 x 10 6 particles). 

In this work, we combine high resolution merger simula- 
tions with exact initial conditions. We present the results 
of a simulation of the massive merger of a 1.1 M and a 



0.9 M carbon-oxygen white dwarf. We follow the evolution 
of the binary system through the merger phase, thermonu- 
clear explosion and nucleosynthesis. Finally we use three- 
dimensional radiative transfer simulations to obtain synthetic 
lightcurves and spectra. 

2. MERGER AND EXPLOSION 

The inspiral and merger i s modeled using a modified ver- 
sion of the GADGET code (ISpringell |2005|). Modifications 
includ e the Helmholtz equation of state (iTimmes & Swestyl 
120001) and a 13 isotope nuclear reaction network that contains 
all a-elements from 4 He to 56 Ni. Radiative cooling effects 
are not included in our simulation. A detailed description of 
the modifications will be given in a forthcoming paper. In 
addition the maximum smoothing length of a particle was re- 
stricted to 10 8 cm. This affects only particles ejected from 
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Fig. 2. — Density and composition of the final ejecta in homologous expan- 
sion 100 s after the explosion of a slice in the x-z plane. The mass fraction is 
shown for carbon, oxygen, silicon, stable iron and radioactive 56 Ni. 

the binary system during the merger but leads to a significant 
speedup of the code. Since these particles are at very low den- 
sities and contain only less than one percent of the total mass, 
they have no noticable influence on the explosion dynamics 
and observables. 

The initial binary system consists of a 1 . 1 M and a 0.9 M 
carbon-oxygen white dwarf constructed from a total of 1.8 x 
10 6 equal-mass particles. Both white dwarfs are set up in iso- 
lation and relaxed with an additional fri ction force for 10 s. 
We then apply the method described in iDan et al.l (1201 ll) to 
slowly move the two white dwarfs close together. When the 
first particle of the less massive white dwarf crosses the in- 
ner Lagrange-point we stop and start the actual simulation. 
At this time, the binary system has an orbital period of about 
35 s. 

The evolution of the binary system is shown in Fig.Q] The 
mass transfer is stable for more than 15 orbits. After about 
600 s the secondary white dwarf becomes dynamically un- 
stable and is disrupted on a timescale of one orbit. As the 
material of the secondary is accreted violently onto the pri- 
mary, material is compressed and heated up on the surface 
of the primary white dwarf. As a consequence hotspots form 
in which carbon burning is ignited. When the first hotspot 
reaches a temperature larger than 2.5 x 10 9 K at a density 
of about 2 x 10 6 s cm 3 we assume that a detonation forms 
(ISeitenzahl et al.N2009h . Note that despite the high resolu- 
tion we use, we still tend to underestimate the temperature 
in the hotspot compared to even higher resolution simula- 
tions (Pa kmor et al.l201 ll) . Only future detailed investigations 
might be able to decide whether or not a detonation really 
forms but the conditions in our SPH simulations suggest that 
it is plausible. 

At this time we map the whole state of the simulation on a 
uniform Cartesian grid with a resolution of 768 3 grid cells and 
a box size of 4 x 10 9 cm. About O.O3M of material is lost by 
the mapping as it is outside the box. Since this material makes 
up less than two percent of the total mass and has a density too 
low to contribute significantly to nuclear burning (i.e. it stays 
unburned) it does not affect the dynamics of the ejecta or the 
synthetic observables derived from the model. 

With these initial conditions we use the LEAFS code 
(iReinecke et al.lfl999ah tha t applies the level- set t e chnique to 
mode l detonation flames (IReinecke et al.l Il999bt iFink et all 
l2010h . We ignite the detonation at the cell with the highest 



TABLE 1 

ANGLE- AVERAGED LIGHT curve parameters for selected bands. 





Bolometric a 
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^max 


18.6 


18.0 


20.8 


23.2 


21.6 


Am\s 


0.74 


1.30 


0.95 


0.67 


0.42 


M(t max ) 


-19.2 


-19.6 


-19.0 


-19.2 


-19.2 


M(W(5)) 


-19.1 


-19.5 


-19.0 


-19.2 


-19.2 



a Note that this is not the true bolometric lightcurve of our calculations. For comparabil- 
ity to observed bolometric lightcurves, the model lightcurve was reduced to UB VRIJHK 
bolometric. 

temperature and follow its propagation through the merged 
object until most of the material is burned. Fig. Q] shows that 
the primary white dwarf and most parts of the secondary are 
already burned 2 s after the detonation formed. The energy 
release from nuclear burni ng unbinds the object. Using an ex- 
panding grid (Ropke 2005)) we follow the dynamic ejecta until 
they reach homologous expansion about 100 s after the deto- 
nation was ignited. They have an asymptotic kinetic energy 
of 1.7 x 10 51 erg. 

3. NUCLEOSYNTHESIS 

In order to obtain detailed isotopic abundances of the ejecta 
we add 10 6 tracer particles to the simulation of the explosion 
that record their local temperature and density. In a post- 
processing step we run a detailed nu clear network contain - 
ing 384 isotopes on these trajectories (Trava glio et al.ll2004l) . 
To mimic the effect of solar metallicity we choose the initial 
composition of the tracer particles as 47.5% 12 C, 50% 16 and 
2.5% 22 Ne by mass. 

In the explosion, a total of 0.7 M of iron group elements 
are synthesized. They consist predominantly of radioactive 
56 Ni (0.61 M ) with a small fraction of stable 58 Ni (0.03 M ) 
and stable 54 Fe (O.O2M ). In addition, 0.5 M of intermedi- 
ate mass elements are produced in the explosion. The ejecta 
contain about 0.5 M of oxygen and about O.15M of un- 
burned carbon. 

The spatial composition and density structure of the ejecta 
in homologous expansion are shown in Fig. O Its inher- 
ently three-dimensional structure is a result of the state of the 
merged object at the time the detonation forms. Since the 
detonation propagates faster at higher densities, the primary 
white dwarf is burned first and its ashes expand. As they ex- 
pand they sweep around the material of the partially intact 
secondary white dwarf that is still being burned. Therefore, 
the ashes of the primary have already expanded significantly 
before burning of the secondary is completed (roughly one 
second later), meaning that the ashes of the secondary com- 
pletely dominate the center of the ejecta. Hence, when the 
ejecta reach the homologous expansion phase several ten sec- 
onds later, the very inner parts of the ejecta do not contain 
material from the primary white dwarf and therefore no iron- 
group elements. 

4. COMPARISON WITH OBSERVATIONS 

Using the three-dimensional density structure from explo- 
sion modeling and the corresponding spatial abundance distri- 
bution from the tracer partic les we apply the Monte Carlo ra- 
diative transfer code ARTIS (iKromer '& Simll2009l: ISimll2007l) 
to calculate synthetic spectra and lightcurves for our model. 

Fig. [3] shows angle-averaged and line-of-sight depen- 
dent lightcurves of our model compared with several well- 
observed normal SNe la. The angle- averaged lightcurves 
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Fig. 3. — Lightcurves of our model. The panels from top left to bottom right contain UBVRIJHK bolometric and broad-band U,B,V,R,I,J,H,K lightcurves. 
The black line corresponds to the angle-average of the model. Grey histograms show lightcurves along seven different lines-of-sight representative for the 
scatter caused by different (100) viewing angles including the most extreme lig htcurves. Th e time is given relati ve to B -band m aximum. The red symbo ls show 
obser vational data of three well-observed normal SNe la, SN 2001el ( Krisciunas et al. 2003), SN 2003du I Stanishev et al. 2007), and SN 2005cf ( Pastorello et al. 

[20071) . 

B-band rise time of 20.8 days, a peak brightness of -19.6, 
-19.0, and -19.2 in the U, B, and V band, respectively, and 
a light curve decline rate in the B band of Ami$ = 0.95). 
These values are w ell within the range of normal SNe la 
dHicken et al.i r2009). More detailed values for selected bands 
of the angle-averaged light curves are given in Table [T] 

In particular, the angle- averaged lightcurves in the U and R 
band agree very well with those observed for normal SNe la. 
The angle-averaged synthetic #-band lightcurve is about 0.3 
magnitudes fainter at maximum than the three supernovae we 
compare with and slightly brighter from 25 days after maxi- 
mum onwards. The angle- averaged V-band lightcurve of our 
model agrees well with the data up to about 30 days after max- 
imum, but declines slightly more slowly afterwards. 

Due to the asymmetric ejecta structure, our model 
lightcurves show a non-negligible sensitivity to line-of- sight 
effects. In the B band, for example, we find peak magni- 
tudes between -19.5 and -18.7. At the same time Am 15 (B) 
varies between 0.5 a nd 1.4, roughly along the Phillips relation 
dPhillips et al.ll 19991) . For the U band, the scatter is even larger 
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Fig. 4. — Maximum light spectrum of our model. The red line shows the 
spectrum of our model one day after maximum light in the B-band. The black 
line shows the observed spectrum of SN 2003du ( Stanishev et al. 2007) at the 
same time. 

show good agreement with observations in the optical bands 
and reproduce general properties of normal SNe la having a 
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while it decreases in the redder bands. For lower photon en- 
ergies, the asymmetry of the ejecta is less important, since the 
optical depths are smaller and photons probe a larger fraction 
of the total ejecta, thereby making the observables less sen- 
sitive to line-of- sig ht effects (compare iKromer & Simll2QQ9l: 
IKromer et al]l2010h . 

The angle-averaged lightcurve of our model in the / band 
agrees well with the observations up to maximum light. After- 
wards it is between 0.2 and 0.7 magnitudes brighter than our 
comparison SNe la. This offset (in the / band after maximum) 
can be attributed to a flux excess in the Ca II NIR triplet in 
the synthetic spectra. This could point to an over- abundance 
of calcium in the model. However, it is more likely to be a 
shortcoming of our radiative transfer treat ment, which uses 
the si mple van Regemorter approximation (Ivan Regemorterl 
119621) to treat collisional excitation. Thus we likely underes- 
timate the effectiveness of cooling by forbidden transitions of 
heavy elements and therefore predict that too much cooling 
occurs via the strong dipole transitions of Ca II. 

There is good agreement between model and data in the 
NIR bands (J,H,K). In particular, our model reproduces time 
and brightness of the second peaks and even the light curve 
decline at more than 40 days after maximum quite accurately. 
However, high precision modeling of the NIR lightcurves re- 
quires an extensive atomic data set to properly simulate flux 
redist ribution by fluorescence ( Kasenl 120061: IKromer & Siml 
120091) . Here, however, we have restricted ourselves for com- 
puta tional reasons to a sim plified atomic data set (cd23_gf- 
5 of IKromer & Siml 120091) with only - 400,000 lines. This 
has been shown to give reliable results in th e optical bands 
but sl ightly underpredicts the NIR brightness (IKromer & Siml 
120091) between the primary and secondary peaks. 

Fig. IH shows the angle- averaged spectrum of our model one 
day after B-band maximum. For comparison we overplot an 
obser ved spectrum of the normal SN 2003du ( Sta nishev et alJ 
120071) at the same epoch. Our model shows most of the char- 
acteristic features of SNe la, particularly the defining Sill 
doublet at AA6347,6371 but also the weaker Sin features at 
AA5958,5979 and AA4128,4131, the Mgii triplet at A4481 
and the Call H and K absorption. The Sll W-feature at 
~5400A (which is a blend of several lines) is also visible 
though it is weaker than in SN 2003 du. In the red tail of 
the spectrum the Ca II NIR triplet is clearly visible as well as 
some indication of the O I triplet AA7772,7774,7775. More- 
over, the overall spectral shape and the velocity- shifts of most 
of the line features in the observed spectrum are well repro- 
duced, indicating that the velocity structure of our model is a 
good representation of that in real SNe la. 

5. DISCUSSION 

The dominant parameter that determines the brightness of 
the explosion in our model is the mass of the primary white 
dwarf. In the merger the primary remains mostly unaffected 
whereas the secondary white dwarf is destroyed. Therefore, 
because the density profile of the primary white dwarf to first 
order only depends on its mass and the density stays high 
enough to burn to iron-group elements ( 56 Ni) only in the pri- 
mary white dwarf, the brightness of the supernova directly 
correlates with the mass of the primary white dwarf. 

Physical parameters of secondary importance are the mass 
of the secondary white dwarf and the composition of both 
white dwarfs (C/O-ratio and metallicity). The material of the 
secondary is burned only to intermediate-mass elements and 
oxygen and therefore only affects the brightness of the ex- 



plosion indirectly. The metallicity of the two white dwarfs 
has only a small impact on the final composition of the ejecta 
by changing the mixture of iron-group elements synthesized 
(ISim et al.ll20lQl) . 

The moment at which the detonation forms is an artificial 
parameter of our model. Physically it will be determined by 
the properties of the binary system discussed above. Since 
we cannot resolve the formation of the detonation microscop- 
ically, however, we have to infer the most likely place and 
time for it to form from macroscopic properties. Once we as- 
sume that a detonation forms on the surface of the primary its 
exact place and time are only secondary effects: the choice 
for the moment of detonation affects the densities at which 
the material of the secondary white dwarf is burned (since the 
detonation occurs while the secondary is being destroyed) but 
it does not strongly affect the burning of the primary because 
it is not dynamically affected by the merger. 

Compared to the merger of two 0.9 M white dwarfs 
(IPakmor et al.ll20ldl) the merger we present here mainly dif- 
fers in the mass of the primary white dwarf. Since it is con- 
siderably more massive, its central density is higher and its 
material is burned predominantly to iron group elements. In 
addition, the mass ratio of the binary system we model here is 
well below unity. Therefore, the more compact primary is not 
subject to noticeable tidal forces from the secondary. 

In addition, although exact initial conditions may weaken 
the merger, we show that it still leads to the formation of 
hot spots on the surface of the primary, in which densities 
and temperatures are sufficiently high that a detonation is 
plausible. This is consistent with the properties of hot spots 
seen in previous studies w ith less relaxed initial conditions 
(IPakmor et al.|[2 010. 20 1 1 ) and also with very recent simula- 
tions by iRaskin et alJ (120111) who found similar hot spots in 
calculations with relaxed initial conditions and on ly slightly 
smalle r resolution than our simulation. Note that iDan et all 
(|2012|) do not find hotspots in similar systems with exact ini- 
tial conditions but more than a factor of ten fewer particles. 

Our model can be interpreted as a pure detonation of an 
isolated sub-Cha ndrasekhar mass white dwarf as described in 
I Sim et alJ (1201 0|) with an additional component from the re- 
mains of the secondary white dwarf surrounding the system. 
Quantitatively, however, the color lightcurves of our model 
agree significant l y bett er with observations than the toy mod- 
els of I Sim et all (120101) . In particular, our model also shows 
good agreement in the NIR-bands and reproduces the position 
of the secondary peaks in these bands very well. 

This difference to the toy models is likely to be associated 
with a different density profile and composition structure in 
the inner parts of the ejecta. As discussed in Section [3] this 
is a result of the presence of the secondary white dwarf and 
its interaction with ashes of the primary when they expand. 
Also, the asymptotic kinetic energy per mass of the ejecta is 
smaller for our merger model than for detonations of isolated 
sub-Chandrasekhar mass white dwarfs, because the burning 
of the secondary white dwarf is less complete and hence the 
average energy release per mass is smaller. 

Overall, our explosion reproduces observational data of 
normal SNe la of the same brightness reasonably well. Since 
the time span from the onset of mass transfer between the 
two white dwarfs until the ejecta reach homologous expan- 
sion is very short (only of the order of a few minutes) there is 
no extended outer envelope the explosion ejecta interact with 
later than a few minu t es afte r the explosion, in contrast to the 
model bv lFrver et all (120101) . 
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Observational constraints dNapiwotzki et al.1 12005L I2007h 
on the number of double white dwarf binaries significantly 
more massive than a Chandrasekhar mass are still not very re- 
strictive. If violent mergers contribute noticeably to the total 
SNe la rate, our calculations suggest that it will be hard to 
distinguish the merger events from the bulk of normal SNe la 
via optical/NIR lightcurves and optical spectra. Further de- 
tailed work is therefore needed to investigate whether there 
are other characteristic properties of SN la explosions from 
violent white dwarf mergers and to explore how the observa- 
tional display changes with the properties of their progenitor 
system. 
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